Born-Infeld corrections to Coulombian interactions 
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Two-dimensional Born-Infeld electrostatic fields behaving as the superposition of two point-like 
charges in the linearized (Maxwellian) limit are worked out by means of a non-holomorphic mapping 
of the complex plane. The changes underwent by the Coulombian interaction between two charges 
in Born-Infeld theory are computed. Remarkably, the force between equal charges goes to zero as 
they approach each other. 



I. INTRODUCTION 



When forces between charges are considered in nonlin- 
ear theories, the picture of the field due to a charge act- 
ing on another charge is no longer applicable. Since the 
superposition principle is not feasible, a multiple charge 
configuration has to be analyzed as a new problem in- 
stead of the mere addition of already known solutions. 
Even the expression "multiple charge configuration" calls 
for an explanation. For those theories behaving as a lin- 
ear theory in the regime of weak field, a multiple charge 
configuration can be defined as a solution of the nonlin- 
ear field equations going to a superposition of individual 
(linear) charges at infinity. Given a static multiple charge 
solution of some nonlinear theory, the force on a charge 
can be worked out by computing the flux of the stress 
tensor through a surface surrounding the charge. Since 
the stress tensor is divergenceless, the force will result 
different from zero only at those points where the stress 
tensor is singular, which provides a way of localizing the 
charges. 

This article is aimed to solve two-dimensional static 
configurations of two charges in Born-Infeld nonlinear 
electrodynamics, and compute the interaction strength. 
In Section II we summarize the Born-Infeld theory. In 
Section III we characterize the two-dimensional electro- 
static solutions by means of a non-holomorphic complex 
transformation. In Section IV we obtain the repulsive 
and attractive interaction between equal and opposite 
charges. We compute corrections to the Coulombian in- 
teraction for distant charges, and show that the repulsive 
force vanishes when equal charges approach each other. 
The conclusions are displayed in Section V. 
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II. BORN-INFELD ELECTRODYNAMICS 



Born-Infeld electrodynamics is a non-linear theory 
whose initial objective was to render finite the self-energy 
of a point-like charge. In Born-Infeld electrostatics the 
electric field E due to a point-like charge does not di- 
verge but goes to a finite value b at the charge position. 
The energy-momentum tensor still diverges at the charge 
position but the the integration of the energy density be- 
comes finite. The fundamental constant b is an upper 
bound for the fields, and regulates the transition to the 
weak field regime: for fields much smaller than b the the- 
ory behaves as Maxwell electromagnetism. By healing 
the field of singularities. Born and Infeld pretended that 
the theory could be regarded from a unitary standpoint: 
the only physical entity would be the field, whereas the 
charges would be just a part of the field P, 0, H, 13, Hi ■ 
They even believed that the solutions will contain some 
essential features of the charge dynamics, which is by 
no means true since the theory allows for static multiple 
charge solutions. 

Born-Infeld electrodynamics possesses outstanding 
physical properties: together with Maxwell theory, they 
are the only spin 1 field theories having causal propaga- 
tion 0, 0] and absence of birefringence d, [1] . Although 
concrete solutions for propagating Born-Infeld electro- 
magnetic waves are not sufficiently known -apart from 
trivial free waves solutions-, solutions for waves propa- 
gating in static background fields and waveguides have 
been recently obtained [§, [l^ . The renewed interest in 
Born-Infeld theory can be traced to its emergence in the 
study of strings and branes: loop calculations for open 
superstrings lead to a Born-Infeld type low energy ac- 
tion [ll|, [12, [3 ■ Nowadays Born-Infeld- like Lagrangians 
have been proposed for quintessential matter models and 
inflation m!4"|, and also for alternative theories of gravity 
[l5| . Born-Infeld charges coupled with gravity have been 
investigated in attempts to remove geometrical singular- 
ities of charged black holes p^ . 

Born-Infeld Lagrangian density for the electromagnetic 
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field Fij = d^Aj - djA, is 0] 
Lbi = (^^\det {bg,,+F,,)\ - ^|det (6.9,,)l 

^b^ (l- Vl + b-^2S~b-^pA (1) 



47r c 



where 5 and P are the scalar and pseudo-scalar invari- 
ants, 



S=\ F,, F^i = H^' - 

E B , 



P = i *Py F'l 



(2) 



and 5 is a new universal constant with units of field, which 
plays the ro le of upper bound for the electrostatic field. 
The term -^/jdet [b gij)\ in ([1]) makes L bi to vanish when 
the electromagnetic field vanishes. Lbi goes to Maxwell 
Lagrangian density in the limit b ^ oo. By defining the 
2-form 



F, , - 



^3 



(3) 



the Euler-Lagrange equations coming from Lbi a-re writ- 
ten as 



(4) 



These equations are supplemented with the identities 
dF = (i.e., Pfjfc.i] = 0), since the field is an exact 
2-form [F = dA) . The energy- momentum tensor is 



which verifies the energy-momentum conservation. 



(5) 



(6) 



at all the places where is non-singular. For a point- 
like charge Q, F = Q r"^ dt A dr diverges at the charge 
position but F is finite. Thus, although the energy den- 
sity still diverges at the charge position, the integrated 
energy turns out to be finite. 

Owing to the non-linear character of the field equa- 
tions ([3]), it is hard to find explicit solutions in an an- 
alytic way, apart from highly symmetric configurations. 
In more general cases, the solutions are displayed only in 
an implicit form [13, [3 • 



III. TWO-DIMENSIONAL ELECTROSTATIC 
SOLUTIONS 

Two-dimensional electrostatic solutions, together with 
conditions guaranteeing uniqueness of solutions, have 
been worked out in a rather cumbersome parametric way 
by resorting to the relationship between minimal sur- 
face equations and the Born-Infeld electrostatic problem 



|19l . |20|. Recently, two-dimensional electrostatic solu- 
tions in Euclidean space have been obtained by using 
a non-holomorphic transformation of the complex plane. 
This method has been used for working out the field lines 
and self-energies of point-like two-dimensional multipoles 
PH . The method searches a coordinate transformation 
in the plane, (x,?/) {u,v), such that (m, u) are orthog- 
onal coordinates and u{x,y) results to be the potential 
for a Born-Infeld electrostatic field: E = — Vu (then the 
coordinate lines are equipotential and field lines respec- 
tively). The equation to be fulfilled by u{x,y) is the one 
resulting from Eq. ^ when the field F is chosen to be the 
exact 2-form F — du /\ dt, corresponding to the electro- 
static potential A = u{x,y)dt. In current language, the 
equation is 



V. 



Vu{x,y) 



= 



\Vu{x,y)Y 



(7) 



which becomes the Laplace equation when b ^ oo. Ac- 
tually this equation can be also obtained by directly sup- 
pressing the time coordinate t and the third spatial co- 
ordinate in the Lagrangian. Thus, in two dimensions the 
(static) field is the exact 1-form F = du. In addition, *JF 
in Eq.© is the 1-form = *P/VT~FW. 

In terms of complex numbers z = x + iy and w — 
u + iv, the coordinate transformation {x,y) — s- {u,v) can 
be regarded as z ^ z(z/;,W) or 



dz = p{w, w) dw + q{w, w) dw 



(8) 



(the bar means complex conjugate). The integrability 
condition for Eq.® requires that 



dp{w,w) = dq{w,w) 



(9) 



where d and d are exterior derivatives with respect to w 
and w (Dolbeault operators). To glance the Euclidean 
metric in {u,v) coordinates, one can write 

dx^ + dy^ =dzdJ^ (IpP + \q\^) \dw\^ + 2 Re{pqdw^) 

(10) 

To get orthogonal coordinates (u,w), the terms contain- 
ing du dv must be ruled out from the quadratic form (fTT 
Then 



Im{pq) = 



(11) 



The 1-form dw can be solved from Eq. ^ and its complex 
conjugate: 



dw — 



pdz — qdz 



(12) 



In two-dimensional Euclidean space ii la * dz ~ i dz\ then 

pdz + qdz 



"dw 



\p\'-k\- 



(13) 



In order that the field = Pe(*dw)/VT~FW ac- 
complishes the Eq.([5]), we have to properly choose the 
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functions p and q. Let us try the choice q = 0. Then, 
the integrabihty condition Q imphes that p — p{w) {w 
is a holomorphic function of z). In this case *dw = idw 
is an exact 1-form. So, Eq.(|3]) is satisfied only if 6 ^ oo 
(because it will result *T — *F ~ Re{*dw)). This is the 
case of the Coulombian field (it is well known that holo- 
morphic functions provide solutions for Laplace equation 
in two dimensions). The Coulombian solution is 



du 
dx 



du 
dy 



1 



p{w) 



(14) 



Let us try a choice of p a nd q leading to g = when 
6 cx). We guess the choice 



p{w)q{w) 



1 
4&2 



(15) 



(notice that the integrabihty and orthonormality condi- 
tions (HI [TT]) are satisfied) . Then 

^ T) I \ . ,2 Re{p) dx + Im{p) dy 

du = Re{w) = 4 h ^^2|p|2 ^ 1 (16) 



dv = Im{w) = 4 -^'^^f^tl^^f^'^'^ (17) 
Therefore 

^^E.^iE.^-^-i^^- j':^^^ . (18) 



dx dy 4 62|p(u;)|2-Kl' 



and 



IF 



± 



4 62 |p(^)|2 



4 62 |p(y,)|2 



1 



(19) 



These results imply that *T is an exact 1-form. In fact 



*T = 



*du 



VI - 6-2£;2 



= ±idv 



(20) 



Therefore, Eq.Q is satisfied. 

Eg. p^ shows the way to generate those Born-Infeld 
field behaving in the weak field limit (far from the 
charges) like a given Coulombian field: choose the holo- 
morphic function p(w) associated to the Coulombian field 
(fH)) and replace it in However, as the field should 

be expressed as a function of z or (x, y) instead of w 
or (m, u), we have to take into account that the relation 
between z and w is no longer the Coulombian relation; 
according to Eq.® the relation now is 



dz = p{w) dw 



1 



Ab^ p{w) 



■ dw 



(21) 



This equation amounts to a non-holomorphic relation be- 
tween z and w. The electrostatic potential fulfilling the 
Born-Infeld equation ([7]) is u{x,y) = Re{w{z,'z)). Func- 
tion p{w) in Ea. (PT|l plays the role of a Coulombian seed 
to obtain the Born-Infeld potential. 



Eg. ITS)) shows that |E| does not diverge but attains its 
maximum value 6 at the points where 



\p{w)\ 



1 

26 



(22) 



According to Eq. (fT5)l . at these points it is \p\ = \q\- so 
they are singular points of the coordinate change p2l [T3| . 
The relation (|22p describes the curve where the energy- 
momentum tensor ([5]) is singular, because the vanish- 
ing of implies that T diverges. This curve is then 
the charge location. If the curve |p(w)| — (26)^^ is 
closed, then it separates two different regions in the com- 
plex plane: i) \p{w)\ > (26)^^ and ii) \p{w)\ < (26)"^ 
Only the first region can realize the Coulombian limit 
E/b 0. Since the Born-Infeld field ^ should go to 
the Coulombian field at infinity, the region \p{w) \ > 
(26)~^ corresponds to the exterior of the charge distribu- 
tion. The curves where ^(u')! = (26) ~^ have been stud- 
ied in Ref. '^21] for the configurations associated with 
Coulombian multipoles. In these cases the curves are 
closed and turn out to be epicycloids whose sizes are de- 
termined by 6 and the multipolar moment. One can say 
that Born-Infeld field smoothes singularities in two differ- 
ent ways: on one hand it smoothes the divergence of the 
energy-momentum tensor in order that the self-energy re- 
sults to be finite; on the other hand the point-like charac- 
ter of the Coulombian multipoles is spread to the surface 
(in this case a curve) where the field reaches the upper 
bound 6. 



IV. FORCE BETWEEN TWO-DIMENSIONAL 
MONOPOLES 

We are going to study the electrostatic Born-Infeld 
configuration corresponding to two point-like monopoles 
separated by a distance d in the Coulombian limit. The 
interaction force between charges will result from the mo- 
mentum flux through a closed surface S containing one 
of the charges. We will choose the x-axis along the line 
joining the charges, and the origin of coordinates at the 
intermediate point. Since the symmetry dictates that the 
force is directed along the x-axis, then the involved mo- 
mentum flux is T^^dSj. As usual, we choose S as the 
surface formed by the y-axis, and a semi-circumference 
of infinite radius centered at the origin of coordinates; on 
this last surface the flux is null. Thus 



= - J T^^ dS 

y — axis 

62 ' 



(23) 



An 



' 62 / 62 V b"^ 



dy. 



x=0 



For repulsive interactions between equal charges it is 
Ex = {v = constant) and dy = (dy/du) du, at x = 0. 
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Thus, according to Eqs. ([THl [HI [H]) the force is 
1 f Im{p{w)) 



Stt J \p{w)\' 

y — axis 



du 



(24) 



For attractive interactions between opposite charges it is 
Ex = ±E {u = constant) and dy ~ (dy/dv) dv, at x = 0. 
Then the force is 



1 



Re ipiw)) 



iTT J \p{w)\' 

y — axis 



dv. 



(25) 



A. Equal charges 

The well known Coulombian potential u{x,y) for the 
repulsive configuration of two equal charges A at a dis- 
tance d can be written as u{x, y) — i?e[w(z)], where w{z) 
is the holomorphic function 



-2A In 



(26) 



In Eq. (pS)) a proper choice of the integration constant was 
done in order that the Coulombian potential be zero at 
the origin. Relation (pS)) can be inverted to obtain the 
holomorphic function p characterizing the Coulombian 
mapping: 



z = ±yexp(-|^ 



1 



(27) 



where ± alludes to x > and a; < regions. By differ- 
entiating the Eq. ([27)) one obtains 



(28) 



d exp 

8A /exp(-i|)+l 



By substituting in Ea. (j2ip and then integrating it, one 
obtains the mapping leading to the Born-Infeld complex 
potential w(z,z): 




(29) 



where a is the non-dimensional parameter a = bd/X. 
In Eq. ((29)) one recognizes the holomorphic (Coulombian) 
seed (|27|) and the anti-holomorphic Born-Infcld correc- 
tion. Differring from the Coulombian case P7|) . Eq.f^H) 
does not provide a unique electrostatic potential u to 



each point of the complex plane. In particular, while 
the Coulombian mapping (|27p is periodic in the v coor- 
dinate, the Born-Infeld mapping ([29]) fails to be periodic 
because of the presence of the linear term W/(4A). Fig- 
ure [1] shows the u, v lines for the Coulombian potential 
([?7)) . The equipotential u ^ passes through the co- 
ordinate origin. The lines v/(2X) = ±7r coincides with 
the y-axis (u > 0) and the part of the x-axis joining the 
charges {u <0). The line = is the piece of the cc-axis 
going from the charges to infinite. Figure [D shows the 
M, V lines resulting from the Born-Infeld mapping (j29p for 
— TT < w/(2A) < TT. Still it is v/{2X) — ±tt at the y-axis; 
however, as a consequence of the lost of the periodicity, 
the lines v/{2X) = ±7r do not end at the x-axis but cross 
it, so giving rise to a multi- valued figure for the potential 
u{x,y). To have a single- valued potential, the u — v do- 
main of the mapping (|29p should be restricted by cutting 
it at the x-axis (owing to this branch cut, the field is not 
continuous along the line joining the charges). Like the 
Coulombian case, the line w = is the piece of the x-axis 
from the charges to infinite. Function p{w) is real on the 
line t; = (see Eq. (j28p '). and attains its maximum value 
1/(2 6) at the charges, where the potential results to be 



(30) 



Since the equipotential lines surround the charges, the 
value ()30|1 is the maximum value for the Born-Infeld po- 
tential. Notice that, by solving the equation |p(w)| = 
1/(2 6) for any value of v one obtains the curve u(v) de- 
scribed by the relation 



-2cos(^)exp(-^) = 



(31) 



1. 



This curve decomposes in two parts at each side of the 
y-axis, displaying cusps on the x-axis at < d/2 (see 
Figure [3|). In spite of the appearances, the charge con- 
figuration is not spread on the lobes of Figure O but 
concentrates at the cusps. In fact, Eq. (|3T|l can only be 
fulfilled for u > Umax- Thus the only admissible solution 
of Eq. (j3ip is {u = Umax, u = 0), i.e. the positions of the 
point-like charges. The rest of the lobes is cut because 
the mapping domain is restricted to have a single-valued 
potential matching the Coulombian potential at the in- 
finity. 

In order to compute the force (124]) we will take into 
account that the points on the y-axis satisfy 



u 

2A < 



2A 



±7r. 



(32) 



Then, according to Eq.f^H), the positive y-semiaxis can 
be parametrized by defining a parameter t such that 
cos(i/2) — exp[u/(4A)]; thus it results 



t A , . ' 

tan — — — (tt + t — sm t) 
2 



(33) 



FIG. 3: Curves |j3(w)| = (26)^^ characterized by a = hd/\ 
for the configuration of two equal charges. Only the cusps, 
but not the lobes, belong to the branch of w{z) under consid- 
eration. 



FIG. 1: Equipotencial and field lines for the Coulombian equal 
charges (right semi-space). V stands for v/{2\). 




FIG. 2: Multi-valued Born-Infeld mapping (I29j for -tt < 
v/{2\) < TT (a^ = 40). In the figure V stands for w/(2 A). 



the Coulombian case (6 oo, i.e. a oo) it is Uo — 0, 
so to = and the Coulomb force is recovered. Instead, 
in Born-Infeld theory it is Uo < and — tt < to < 0. The 
dependence of the interaction on h is given through the 
value to which depends on a = bd/\. Since to — sin to is 
negative for — tt < to < 0, it is concluded that the inter- 
action between equal monopoles is repulsive for all value 
of a, but less intense than the Coulombian interaction. 
In particular, for a — > 0, it is to — tt; thus the repulsive 
interaction force vanishes for d 0. 

In order to get an explicit correction to the Coulombian 
force we will take into account that to is small for high 
values of a. Thus we can try to solve y{to) = in ([33)l by 
writing to as a power series of a~^. In this way we reach 
the result t, 
O(a-iO) 



Therefore the force ([M)) is 
256 TT^^ 



2A2 



1 



3 



+ 0{a-'^). (35) 



Then, Born-Infeld correction to the repulsive force be- 
tween equal charges is very weak. Notice that d in ([55)1 
is not the real distance D between the charges. D goes to 
d when 6 ^ 00, but D is smaller than d for equal charges 
(see Figure |3]): 



D 



repulsive ~ ^ ^"^i.^) ^ l) 



0\a 



(36) 



where t e (— 7r,to). Then the integral (I24p can be com- 
puted as twice the integral along the positive y semi-axis. 
Thus the force is 



2A . 


^ exp 


" U ' 




U ' 


~d J 




.2A. 


y^exp 


.~2A. 



1 du 



—00 

to 



= — / sm^ T;dt = —7 (tt + to - sm to) (34) 
TT a J 2 TT a 

— 7r 

In the former expression, Uo is such that y{uo) = 0. In 



B. Opposite charges 

We will now repeat the former steps for the case of the 
attractive charge configuration consisting of a charge — A 
at a; = — d/2, and a charge A at a; = (i/2. The complex 
Coulombian potential is 



w ■ 



-2A In 



2z _ 1 
d ^ 



2z_ 

d 



1 



(37) 



Here we have chosen the integration constant such that 
w{z) is null at infinity. Inverting (j37p we obtain the 
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Coulombian mapping 

d w 
z - 2 4A ' 

This means that function p{'w) is 

/ \ ^ ■ 1^-2 w 

P(w) = smn — 

' 8A 4A 

Thus the Born-Infeld mapping results 
d 



'=2 



cotn — H T tt; — smli — 



4A 



2A 



2A 



(38) 



(39) 



(40) 



The points where the electrostatic field reaches the value 
b (i.e., |p(it;)| — (2 6)^^) belong to the curve 



u 
2A 



± arccosh 



(2+^°^ 2aJ 



(41) 



If a > 4 this curve decomposes in two separate parts 
("charges") at each side of the y-axis. Like the previous 
case, the Born-Infeld mapping fails to be periodic 
due to the presence of a linear term; so the u — v do- 
main in mapping (I40p should be properly restricted to 
get a single-valued potential u{x,y). The branch of the 
complex potential w to be kept is those matching the 
Coulombian potential at infinity. Differing from the pre- 
vious case, this branch does reach the curves (|4ip . Fig- 
ures m and [5] show the equipotential {u = constant) and 
field {v = constant) lines surrounding the right charge, 
and their relation with curve (|4T|) . Again the lines v = 
coincide with the piece of the x-axis going between the 
charges and infinity. But in this case the domain of 
|t;|/(2 A) has to be extended beyond tt to reach the piece 
of the X-axis joining the charges. Both the potential u 
and the field E are discontinuous at the charge (the field 
attains the maximum value b at the exterior side of the 
charge, i.e. the side where < |'y|/(2A) < tt). Besides, 
the field is discontinuous at the branch cut on the x-axis 
between the charges. 

For a < 4 the curve of maximum field get closed, as 
it is typical for a dipole 2l|. In this case, the "charge" 
distribution becomes a unique object; the Coulombian 
zone |p(w)| > (26)~^ is, of course, the outside of the 
object. Figure [S] shows the curves |p(w)| = (26)"^ for 
different values of a. The curves display cusps on the 
X-axis {v — 0) at |a::| > d/2. At the cups it is u = 
±2A arccosh (q;/2 + 1). 

We will calculate the force (|25p for a > 4, which is 
the case where the charges are separate. Since the j/-axis 
is characterized by m = 0, then Eq.fgDl) says that the 
positive y-axis accepts a parametrization similar to (I33p 
whenever the parameter t is defined as t = t;/(2A) — tt: 



yit) 



tan ■ 



t 



;{n + t + sini) 



(42) 



This function is monotonous for t € (— tt, i^), where 
< to < TT is the parameter satisfying y{to) — 0. Since 




FIG. 4: Born-Infeld equipotencial lines (u — constant) for 
a — 5, together with the curve describing the right charge. 
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— v=^i'^^x\ \ 
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O.iS 0.5 


0.75 
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l.i5 1.5 


1.75 
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FIG. 5: Born-Infeld field lines [v — constant) for a — 5, 
together with the curve describing the right charge. In the 
figure V stands for v/{2X). 



)\y-ax^s ^ d/{8X) siu"^ b/(4A)], then the force ^ 

is 

to 

4 A^ /" t 2 A^ 

= / cos^ -dt = {Tr + to + sin to) (43) 

wd J 2 nd 



For 5 — > 00 (a ^ 00), it is to ^ in Eg. (1^5]) . so the 
Coulombian force is recovered. Since to € (0, tt), it is con- 
cluded that the attraction between opposite monopoles 
is more intense than Coulombian interaction. We will 
solve y{to) = in by writing to as a power series 
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y/d 



^P;. 6^-0.4 -0.2 

\ ^ - ^ 



a= 5 - - ^ \ \ 



x/d 



0.2 0.4 



FIG. 6: Born-Infeld opposite "charges" for a = 10 and a = 5. 
For a < 4 the charges merge in a sole object. 



of (o? - 16)"^ The result is to = Stt {c? - 16)"^ (1 - 
167r2 (a^ - 16)^73) + 0[(a2 - 16)"'^] . Therefore the Born- 
Infeld interaction (^5]) between opposite monopoles be- 
haves as 



2A^ 



16 



512 TT^ 



16 3(a2 



- 16)3 
16)-* 



(44) 



Differing from the repulsive case, the attractive inter- 
action receives a more perceptible correction of order 
(o? — 16)-^. Notice that d in is not the distance 
D between the cusps in Figure [SI By computing the po- 
sitions of the cusps for opposite charges, it results that 
D is larger than d: 



D. 



attractive 



= d 



2d 



or — 16 



(ln(a^ - 16) 



+ 0[(a^-16)-3/2] . (45) 



V. CONCLUSIONS 

The first goal of Born-Infeld theory was the obtain- 
ing of a point-like charge solution with finite self-energy. 
However this solution has intriguing features: dLsi/d'E 
still diverges, and E is finite at the charge position (an 
unpleasant property for a vector field at its center of 
symmetry). However, these disagreeable features do not 
cause any trouble to the interaction between charges. We 
have considered a "two charges" field (which differs from 
the mere superposition of two "one charge" fields, since 
the theory is nonlinear). To compute the interaction be- 
tween parts of this field configuration one must consider 
the momentum flux through a surface separating both 
subsystems. According to the method developed in Ref. 
[2l| . the interaction force is given by expressions (17- 
18), where p is a Coulombian function, and the Born- 
Infeld features are encoded in the integration interval. 
Although we have obtained the force in a parametric form 
(the parameter to in forces (I34p and (I43p comes from the 
transcendent equation y{to) = in ([33l) and ((42|) respec- 
tively), we have succeeded in computing Born-Infeld cor- 
rections to Coulombian interactions. In addition we have 
proved that the interaction force between equal charges is 
well behaved and goes to zero when the charges approach. 
This limit cannot be reached for opposite charges because 
they merge in a unique dipolar object of finite size. 
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